% 测试1：生成标准正弦波（幅值1，频率10Hz，采样率100Hz）
fs = 100;
f0 = 10;
t = 0:1/fs:1;
x = sin(2*pi*f0*t);
Q = fourier(x, fs, f0); % 期望值 Q ≈ 0.5（因为正弦波FFT幅值为N/2，Q=(N/2)/N=0.5）
disp(Q);

% 测试2：超频警告
try
    Q = fourier(x, fs, 60); % fs=100Hz时，Nyquist频率为50Hz，60Hz超限
catch ME
    disp(ME.message); % 应触发错误提示
end